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ABSTRACT 

We analyse behavior of stellar winds of evolved stars and the outcome of SN la explosions 
in the AGN driven relativistic bubble. We find that the expansion of wind shells is efficiently 
decelerated by the relativistic pressure; their bulk motion however is preserved so they cross 
the bubble together with the parent star. The wind material occupies a small fraction of bubble 
volume and does not affect substantially the expansion of SN remnants. The estimated maxi- 
mal radius of a SN remnant in the bubble is 30-40 pc, if the envelope keeps its integrity and 
remains spherical. A fragmentation of SN shell due to Rayleigh-Taylor instability can allevi- 
ate the propagation of the SN material so the ejecta fragments are able to cross the relativistic 
bubble. Outside the bubble wind shells and supernova fragments are decelerated in the intr- 
acluster medium at close range off the bubble boundary. The deposited SNe la material can 
enrich the intracluster gas with metals in a thin layer at the boundary of the relativistic bubble. 
This process may lead to a rim of enhanced line emission. In the opposite limit, when the 
fragmentation of supernova remnant is moderate or absent, the SN la matter is advected by 
the relativistic plasma and may leave the central region of the bright cluster galaxy together 
with buoyantly moving bubbles. 
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1 INTRODUCTION 

Radio and X-ray observations show that in majority of regular 
galaxy clusters, possessing a cool core, the activity of a cen- 
tral supermassive black hole mediated by AGN jets creates bub- 
bles of relativis t ic pla s ma in the intra-cluster medium (ICM) 
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consistent with the assumption that bub bles are completely devoid 
of thermal gas jSanders & Fabianll2007l) although the limits on the 
amount of thermal gas in the bubbles are not very tight. The ab- 
sence of strong shocks in the ICM surrounding the bubbles implies 
that the relativistic plasma is in approximate pressure equilibrium 
with the ICM. These bubbles, inflated by an AGN, are believed to 
be responsible for mechanical coupling of the AGN energy release 
and the thermal state of the ICM in galaxy clusters, groups and in- 
dividual elliptical galaxies. For this reason every aspect of bubble 
physics receives much attention. 

Given the ubiquity of bubbles in clusters, these objects should 
be almost always present in the cluster core. The lifetime of the 
bubble in the stratified atmosphere is set by "buoyancy" time, while 



the g rowth rate of the bubble is defined by the pow er of the AGN 
(e.g.. lGull & NorthovedfT973l : IChurazov et al.l l2000l). Pairs of bub- 
bles are located on both sides of the central supermassive black hole 
and often several generation of bubbles differing in their size and 
the distance from the center are observed. 

Bubbles with sizes of order 1-10 kpc are often found in the 
inner regions of galaxy clusters, sharing the space with brightest 
cluster galaxy (BCG). For example, in NGC 1275 (BCG of Perseus 
cluster) within the effective radius of the galaxy ~15 kpc) we 
see two bubbles on either side of t he nucleus with the ra dius of 
each bubble ~ 6.3 kpc (e.g., Bohriii ger et al.i ri993: Fabi an et"al] 
'200: ^. In M 87 (Virgo clii ster) the bubble radius is ~ 1.4 kpc (e.g., 
iForman et al.ll2005l |2007|) , while the galaxy effective radius is of 
order 7.7 kpc. We expect that in both cases bubbles should contain 
significant fraction of the galaxy stars. These old low mass stars 
evolve as usual, lose their mass via the wind, and some of them give 
birth to type la supernovae (SNe la). Even if the relativistic bubble 
sweeps up all the thermal gas in the process of inflation a lot of 
material in the form of the stellar wind thus can be supplied by the 
stars inside the bubble during the bubble life. This wind material 
could affect the expansion dynamics of supernovae exploded in the 
bubble. 
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The question we address here is what happens to the stellar 
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winds and SNIa ejecta embedded in tlie relativistic bubble. Almost 
weightless relativistic plasma provides highly unusual conditions 
for the dynamic evolution of wind shells and SNe. Two extreme 
scenarios are conceivable. In one limit all the wind material and 
SN ejecta are decelerated and mixed with relativistic plasma. Dur- 
ing subsequent evolution the buoyantly rising bubbles advect this 
material with them. In another limit the matter ejected by evolv- 
ing stars propagates freely through the relativistic plasma, attains 
the boundaries of the bubble and enriches the ICM just outside the 
bubble with heavy elements. We investigate different scenarios of 
dynamical evolution of wind shells and SN ejecta and explore ob- 
servational outcomes of these scenarios. 

The structure of the paper is as follows. We first study the 
expansion and bulk motion of the wind envelopes including mass 
striping efi^ects. This provides us with the estimate of the filling 
factor of the wind shells ensemble. We then address the issue of 
the SN envelope expansion dynamics and bulk motion, Rayleigh- 
Taylor fragmentation of the SN shell and propagation of ejecta frag- 
ments in the relativistic plasma. Finally, we discuss the implications 
of results for the relativistic bubble matter contents and intracluster 
thermal environment. 



2 WIND MATERIAL IN RELATIVISTIC BUBBLE 
2.1 Stars and mass injection rate 

We consider, as a fiducial model, a spherical bubble of radius Ri, = 
5 kpc with the bubble center at the distance Ri, from the center of the 
BCG. The characteristic age of this bubble is ti, ~ IRiJv), ~ 3.3 x 
10^ yr, where vve adop t the bubble rise velocity vj, = 300 km s"' (cf. 
IChurazov et alj200lh . At the radii r < 15 kpc t he stellar componen t 
in massive elliptical galaxies dominates (e.g.' Johnson 
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and its density distribution can be approximated by the singular 
isothermal sphere with the velocity dispersion cr,, 
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where G is gravitational constant and ro is a radial scalfl Adopt- 
ing rp = 10 kpc an d velocity dispersion tr, = 350 km s"' 
dWu & Tremainell2006h one gets the stellar mass M{< ro) = Mq = 
5.6 X lO" Mq. This value is consistent with estimates of the to- 
tal mass of stellar component of BCG of ~ 10'" Mq. Inside the 
bubble of radius Ri, = 5 kpc the stellar mass in this case is 
a 9 X W^°R,,_s Mo, where R^.s = Rh/(5 kpc). 
The stars are presumably old with an age oft - 10'° yr, which 
suggests that th e current upper limit of the stellar mass at the AGB 
stage is « 1 Mo ( ISchalleretalJl992l) . Assuming the Salpeter initial 
mass function dN/dm = CnT" within the range mi < m < nii one 
gets the normalizing factor 



C = (a-2)M,[m. 



: 0.28M, , 



(2) 



where a = 2.35, mi =0.1 Mq, and m2 = 1 Mq are used. 

Th e star of m = 1 M p leaves behind a white dwarf of m,,,^ = 
0.5 Mo jSalaris al.ll2009l) . while Am - m - m„ii = 0.5 Mq is lost in 
the form of slow (u x 10-30 km s~') wind during the thermally 
pulsing AGB stage t Vassiliadis & WooJl993h . The present day in- 
tegrated rate of the wind matter injection into the relativistic bubble 
is 



' This approximation breaks at large radii to ensure convergence of the 
total stellar mass. 



where dnij/dt is the rate at which upper limit m2 decreases with 
time. This rate is determined by the relation between the lifetime 
and initial st ar mass t = tiimiln if, where p = 3.2 in the range 
of 1 - 2 Mo JSchalleret aLlll992h . For these values and M, = 9 x 
10"'/?fc,5 Mo the equation 0} yields M = 0.416 M^Rt^s yr"'. The 
corresponding stellar death rate is A' = O.SSRh.s yr"' which is also 
the rate of wind shell formation A'^. One expect thus to find N„, = 
tbN„ X 2.7 X lO' Ri,^ newly created wind envelopes in the bubble 
volume with the total amount of the wind matter in the bubble of 
0.5A',, a 1.4 X 10'/?i.5 Mo, factor ~ lO* larger than the mass in 
the form of relativistic particles of the bubble, 3pV/c- ~ 2.4 x 
10'' Mo. It should be noted, however, that the A',,, estimate ignores 
so far a possible escape of wind envelopes from the bubble, which 
is addressed below. 



2.2 Wind shell dynamics 

The dynamical elfect of the wind matter in the bubble on a certain 
wind shell or SN ejecta depends on the filling and covering factors 
of wind shells. To assess the situation one needs first to find the 
average volume and size of a wind envelope at the final stage of its 
expansion. Here we assume that the bubble is static and adopt that a 
star moves with the mean velocity v,. Despite a singular isothermal 
sphere is assumed for the stellar population, it is reasonable to esti- 
mate Vj using Maxwell velocity distribution truncated at the escape 
velocity v^. T he truncated Maxwell distribution is taken in the form 
proposed by ( lKing|ll966h : /(v) oc [exp(-v'/2o-;) - exp(-v2/2o-;)] 
for V < Vg and / = otherwise. Adopting cr,, = 350 km s"', i.e., 
= 2cr = 700 km s"' we come to the average velocity v, = 400 
km s"'. 

The geometry of the wind shell that forms as a result of the 
interaction of the wind with the relativistic medium depends on the 
value of the drag force exerted on the wind boundary. If the drag 
force is strong the stripped wind material creates a trailing plume. 
On the other hand, if drag force is very weak, then the wind shell 
moves with the star velocity retaining spherically-symmetric shape 
and eventually may escape the relativistic bubble. To explore this 
issue we conside r major stages of the mass loss at the AGB and 
post-AGB stage: jSteffen et al.lll998l : lLou et alj20I(]t) : (1) the slow 
wind (u = 10 km s"', M ~ 10"^ Mo yr"') on the time scale of the 
AGB stage, i.e., 10* yr, (2) slow superwind (M ~ 10"^ - 10"" Mo 
yr"') during the last ~ 10"* yr of the AGB stage, and (3) fast wind at 
the post-AGB stage which corresponds to the planetary nebula (PN) 
stage (~ 10* yr). The last stage practically does not contribute to the 
mass loss, but turns out essential for the acceleration of the slow 
wind. We adopt that 60% of the hydro gen shell is lost at the first 
stage and 40% at the superwind stage jSteff'en et al.lll998l) . which 
implies the mass-loss rates M ~ 3 x 10"' Mo yr"' at the slow wind 
stage and M ~ 2 x 10"^ Mo yr"' at the superwind stage. For the 
fast wind stage we adopt parameters derived from the modelling of 
the X-r ay emission of th e PN: M ~ 2 x 10"* Mo yr"' and u = 1500 
km s"' jLou et al.l201(]h . The fast wind parameters suggest that the 
total kinetic energy released during this stage (~ 10* yr) is £3 » 
4.5 X 10*^ erg. When transfered to the slow wind shell with the mass 
of 0.5 Mo this energy accelerates the shell up to 11 x 30 k m s"', in 
accord with the expansion velocities of evolved PN (Ric her et al.l 
I2OO8I) . 

The mass striping rate of the wind shell scales as the shell 
radius squared, so at the stage of slow wind the maximal stripping 
is attained at the end of this stage. The radius of the wind shell at 
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this age can be estimated from the energy arguments. The kinetic 
energy of the wind shell is spent on the pV work against external 
pressure p and on the internal energy which results in the stopping 
radius of the wind shell 

1/3 
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where pio = p/ClO"'" erg cm"^), y = 5/3, ui = 10 km s"', and 
Ml = 0.3 Mq are used. 

The wind shell moves as a whole together with the white dwarf 
unless the bulk of the material is stripped into the trailing plume. 
The wind shell stripping can be estimated following Nulsen ( 1982) 
consideration of the gas stipping for a galaxy moving in the ICM. 
The turbulent stripping is determined by the combined effect of 
the Kelvin-Helmholtz instability (KHI) and the ram pressure drag. 
Indeed, KHI broadens the boundary layer which results in the su- 
pression of the KHI. The net stripping rate, therefore, is controlled 
by the ram pressure 
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where we use pa = 3 pic- for the density of ambient medium. The 
average residence time of a wind shell in the bubble is i?;,/v, ~ 
1.2 X 10' yr, so the above mass loss rate implies that the wind shell 
loses ~ 2 X 10"' Mq while moving in the bubble, negligibly small 
amount compared to the mass of the wind shell, 0.3 Mq. 

Alternatively, the stripping could be caused by the Alfven 
wave drag. In this regard we note that the infinite conductivity ap- 
proximation for the wind shell is fully applicable. A conducting 
body moving with velocity v across the magnetic fiel d B experi- 
ences the drag force due to Alfven wave generation jPrell et al.l 
[l96l 
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where 5 is the area of lateral surface perpendicular to S, is the 
Alfven velocity, a; S/ ^jAixpa with pa = 3p/c-. Strictly speak- 
ing, the Alfven velocity in the relativisic plasma ( lGedaliiilll993l) 
is smaller compared to this expression by a factor of 0.7-0.9 de- 
pending on the ratio of magnetic to total pressure; we neglect this 
difference. For a sphere of the radius r the lateral area is 5 » 2;rr' 
and the mass stripping rate caused by the Alfven wave drag is 

d2,.2 

O . - - ~-10 7/6, 



M= — ^ «4.1x 10-'V;^°S5/Woyr-', (7) 
Zva 

where = S/(10"' G). The stripping rate due to the Alfven drag 
thus turns out two orders of magnitude larger than the rate accord- 
ing to equation l|5). Yet even for the Alfven drag the mass lost dur- 
ing the residence time is only ~ 0.01 Mq which is a small fraction 
(~ 3%) of the wind shell. We thus conclude that the wind shell 
lost at the slow wind stage remains almost intact while traveling 
across the bubble. A similar result can be obtained for the super- 
wind stage. The outcome of a combined effect of all three stages of 
the mass loss is a spherical wind shell of 0.5 Mq expanding with 
the velocity of 30 km s"'. Using equation l|4j one finds the wind 
shell stopping radius is r„, = O.Sp'^g^ pc. 

The above treatment of the wind dynamics suggests that the 
cosmic rays diffusion into the shell can be neglected. To check 
whether this assumption is valid we assume Bohm diffusion coef- 
ficient rgC/3, where rg is the proton giroradius. This assumption is 
standard for the analysis of cosmic ray propagation and supported 
by observ ational data on th e cosmic ray acceleration in supernova 
remnants jStage et alj2006h . although the concept of a tangled field 
might seriously modify a picture of the cosmic ray diffusion in 



magnetic field jNaravan & Medvedevll200l]) . For the spectral in- 
dex of relativistic protons > 2 the energy of cosmic rays resides 
in low energy protons, which permits us to use the characteristic 
energy of relativistic protons ~ 1 GeV. With B = 10"' G one gets 
r,, ~ 3 X 10" cm. The diffusion time is then 
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where r„, ig = r„,/(10'** cm). For the slow wind stage with the final 
radius of ri ~ 0.2 pc the diffusion time is comparable to the dura- 
tion of this stage (~ 10'' yr), while the life time of superwind and 
fast wind stages is significantly smaller than the diffusion time. 

The estimated time scale of the cosmic ray diffusion suggests 
that the penetration of cosmic rays in the wind can affect the wind 
expansion dynamics at the slow wind stage making the pressure 
gradient smoother and the deceleration less pronounced. As a re- 
sult, the final radius of the wind shell at the slow wind stage in fact 
could be somewhat larger, ri > 0.2 pc. On the other hand, the effect 
cannot be significant because the diffution time increases oc so 
the role of the cosmic ray diffusion rapidly drops for larger radius. 
We conclude therefore that the stopping radius of the wind shell 
(r„, ~ 0.5 pc), which includes the combined effect of AGB and post 
AGB mass loss and omits the cosmic ray diffusion, is a reasonable 
estimate. 



2.3 Wind shell escape 

Outside the relativistic bubble the wind shell turns out in the in- 
tracluster thermal gas. For the mass stripping rate M = nr-f^paVs 
adopting the wind shell radius r„, = 1.5 x 10'^ cm, number density 
of interstellar gas n = 0.02 cm"', and v, = 400 km s"' one obtains 
M ~ 1.7 X 10"' Mq yr"'. It takes 3 x 10* yr to completely strip the 
0.5 Mo wind shell over the distance of ~ I kpc. The same estimate 
can be obtained from the equation of deceleration of the wind shell 
as a whole by drag force ;rr^,p„Vj. The wind shell escaping bubble 
is decelerated thus in a close vicinity of the bubble boundary. 

The average residence time of the wind shell in the bubble 
Rb/^s ~ l -2x 10' yr is somewhat smaller than the age of the fiducial 
bubble 3 X 10' yr. The total amount of the wind shells residing in 
the bubble is, therefore, NRjjjv^ ~ 10', while the filling factor of 
the ensemble of wind shells in the bubble is 



/ = KXRJRt? ~ 10-' p]^ . 
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The probability of a collision with the wind shell is determined by 
the ratio of the bubble radius and the mean free path. The latter is 

A = (7Trln,,y^ = 64pj(^' kpc 

where n„. = (3/47t)N„./RI = 2 x 10"' pc 
of wind shells. The probability of shell collisions is low, because 
the average number of wind shells along the bubble radius is only 
T = Rh/A = 0.08. The average probability of the collision is ap- 
proximately a [1 - exp(-T)] » T = 0.08. More accurate estimate 
can be done using expression for the escape probab ility of the pho- 



(10) 

' is the number density 



can be done using expression tor tne escape probab il 
ton from a homogeneous sphere ( IOsterbrocM[l989l) 



(11) 



For T = 0.08 the equation Jl lb gives pesc = 0.94. Most of wind 
shells therefore escape the bubble freely and only 6% of wind shells 
collide with another shell. 

The fate of the collided wind shells depends on whether the 
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Figure 1. Thin shell model of supernova evolution in relativistic bubble. 
Shown are the shell radius (a), evolution of velocities of the shell (S), pre- 
shock velocity of supernova ejecta (SN), and reverse shock speed (RS) (b), 
evolution of the reverse shock temperature (c), X-ray luminosity (d). 



collision is adiabatic or radiative. With the average relative ve- 
locity of collision u ~ 560 km s"' and the wind shell density 
~ 200 cm"' the estimated cooling time of the shocked gas tumes 
out to be ~ 3 X 10'° s, comparable with the hydrodynamic scale 
r„lu ~ 3 X 10'° s. This means that both adiabatic and radiative 
collision regimes are plausible. In the adiabatic case wind shells 
approximately retain their sizes and absolute velocities, so the adia- 
batic collision does not affect their escape. In radiative case collided 
shells merge and form thin dense pancake of a thickness « r„ and 
density ~ {r„lb)p„, where p„, is the density of the wind shell be- 
fore collision, This pancake is liable to fragmentation into clumps 
of size a > b and average velocity ~ v,/ V2. It is easy to verify 
that for fragments with sizes a > b and density ~ (r^r/h)p„, the 
stripping (or deceleration) time is greater than the stripping time of 
wind shells. We thus conclude that the most of the wind material 
escapes into the hot ICM. 

In our analysis of the wind shell dynamics we ignored a possi- 
ble fragmentation of the wind shell due to the Rayleigh-Taylor (RT) 
instability on the deceleration or acceleration stages. This omission 
facilitates the consideration; yet it does not affect the major conclu- 
sion that the wind shell material escape the relativistic bubble. As 
we will see below, the RT fragmentation favours the shell matter 
escape. 



3 TYPE lA SUPERNOVAE IN RELATIVISTIC BUBBLE 

With SN la production efficiency if/ = 0.008 per one white dwarf 
formed in the stellar population of E-galaxies jPritchet et al.ll200^ 
and the stellar death rate in the bubble of fiducial model A' = 0.83 
yr"' one expects if/Nti, ~ 2 x lO"" SNIa explosions during the rela- 
tivistic bubble lifetime = 3 x 10' yr We now consider in detail 
SN expansion in the relativistic bubble and analyse an outcome of 
the Rayleigh-Taylor fragmentation of decelerating SN shell. 



3.1 SN expansion 

Given a small number of wind shells along the bubble radius (Sec- 
tion |2.3l l the deceleration of the SN expansion in the relativis- 
tic bubble is probably dominated by the pressure of the relativis- 



tic fluid. Indeed, for SN expanding with the characteristic veloc- 
ity V a (2£'/M)''- = 10' cm s"' one readily sees that the ram 
pressure is small compared to the pressure of relativistic fluid: 
pv" = 3p(v/c)^ « p. The crucial role of the external relativistic 
pressure in the SN deceleration is a distinguishing feature com- 
pared to the standard case of SN shell in the ordinary interstellar 
medium. 

In our analysis of the SN expansion we assume the isotropic 
pressure of the relativistic medium. This is the case, if the mean 
free path for relativistic protons along the magnetic field is much 
less than the SN radius. Since the SN expands subsonically relative 
to the external medium, in which the sound speed is « c/ Vs, the 
strong forward shock does not form. The reverse shock obviously 
forms, because outer layers of ejecta are decelerated by the external 
pressure and the velocity jump between the undisturbed ejecta and 
swept-up shell exceeds the sound speed in the unshocked ejecta 
(~ 10 km s"'). The SN is fully decelerated when the reverse shock 
crosses the bulk of the ejecta mass. 

To estimate the stopping radius of SN one can use the energy 
considerations likewise we did for the wind shell expansion. The 
initial kinetic energy E of SN should be spent on the PV work 
against the external pressure and on the internal energy pV/(y - I) 
which gives the stopping radius 
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With a characteristic ejecta velocity v » 10' km s"' it takes fj ~ 
Tin/v ~ 3 X 10' yr to reach ?■„,, rather short time compared to the 
bubble lifetime. 

The dynamics of the swept-up shell and the X-ray emission of 
the reverse shock can be illustrated using a model based on a thin 
shell approximation. This suggests that the shell formed by ejecta 
material flowing into the reverse shock is considered as a thin shell 
which dynamics is governed by the dynamical pressure of the SN 



ejecta and external pressure p 
motion for the thin shell is 



10 erg cm . The equation of 



dv , 
M — = 4nr 
dt 



where the shell mass is determined by the mass conservation 
dM 
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dt / 
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These equations are solved numerically assuming a freely expand- 
ing SN with the mass of 1.4 Mq, energy of 1.5x 10^' erg, the density 
distribution p oc exp (-v/vo), boundary velocity of 4 x 10'* km s"', 
and initial outer radius of 10'** cm. 

Results are displayed in Fig. [T] which shows the evolution of 
the shell radius, velocity of the shell, boundary velocity of SN 
ejecta and velocity of the reverse shock, reverse shock temperature 
assuming full equilibration, and X-ray luminosity of the reverse 
shock. The maximal radius of the thin shell model is 39 pc, slightly 
larger than analytical estimate r;„ = 36 pc. Remarkably, the thin 
shell shows contraction phase (Fig. la) which is a direct outcome 
of the dynamical role of the external pressure. However, since we 
neglect the internal pressure of the shocked envelope, the amplitude 
of the contraction phase in our model is exaggerated, so we stop the 
computations at this phase. To calculate the X-ray emission we as- 
sume that the hot plasma in the shell is distributed homogeneously 
in the shell with the thickness AR/R = 0.1. It is rather a crude ap- 
proximation because the density distribution in the reverse shock is 
expected to be essentially inhomogeneous with a peak at the con- 
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Figure 2. Density distribution for SNIa envelope expanding in tlie ICM at 
the ages -60, 400, 3000, 2 lO'' and 10^ years for tlie ICM temperature of 
0.01 keV (bottom) and of 100 keV (top). Position of a sharp wiggle in the 
density distribution con'esponds to the contact discontinuity separating SN 
ejecta and the ICM. 
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Figure 3. Radius of the contact discontinuity, separating SN ejecta and the 
ICM as a function of time. Four curves shown correspond to explosions in 
the ICM with the same pressure, but different temperatures, 0.01, 1, 100, 
10* keV, from bottom to top. Clearly the final size of the ejecta is largest in 
the ICM with highest temperature. 



tact surface. Yet it is reasonable enough to get an idea about X-ray 
luminosity within a factor of two. The luminosity is maximal at 
the contraction phase and reaches ~ 4 x 10^^ erg s"' at the shock 
temperature of ~ 200 keV. The equilibration of electrons and ions, 
however, is an oversimplification, so the shock electron tempera- 
ture and the luminosity should be considered approximate. 

Another interesting view on the SN expansion dynamics in the 
relativistic bubble gives us one-dimensional hydrodynamic simula- 
tions in which we assume a hot rarefied non-relativistic plasma to 
be a proxy for the relativistic fluid. The thermal pressure is taken 
the same p = 10"'" erg cm"\ The ICM temperature varies in differ- 
ent runs from 10"^ up to 10"* keV0. For 7" = lO'' keV the situation 
is close to the case of the relativistic medium because the thermal 
pressure exceeds the dynamical pressure in the upstream flow of 
the forward shock. We assume homologous expansion of the enve- 
lope V cc r and model the initial density distribution 10 years after 
the explosion as p oc e"'''", where ro = 3 10"^ pc. The ejecta mass 
is 1.4 Mq and kinetic energy is 1.5 x 10^' erg, while the maximum 
expansion velocity is set to 2 10* km s"'. 

The dependence of the expansion dynamics on the tempera- 
ture of the ICM (at the same pressure) is apparent from Fig. |2l In 
the low temperature case (bottom panel in Fig. |2j most of ejecta 
energy is spent on a forward shock, which is barely resolved in our 
simple simulations. By contrast, for the high temperature ICM al- 
most all the initial kinetic energy is eventually converted into the 
enthalpy of the ejecta and only tiny amount of energy is deposited 
in the forward shock. Accordingly the final size of the envelope at 
the boundary separating ejecta and ICM is much larger for the high 
temperature run. This is further illustrated in Fig. [3] showing the 
time dependence of the envelope radius. The simulations show that 
in the limit of very hot ICM the final envelope radius converges 



We use nonrelativistic equation of state in these illustrative runs even for 
Tf = 10* keV, although this not valid for electrons. 



to the value given by the thin shell model (Fig.[Tll. The X-ray lu- 
minosity of the reverse shock in the model with 10* keV medium 
is shown in Fig. |4] together with the evolution of the radii of the 
reverse shock and contact surface. The luminosity evolution is con- 
sistent with the prediction of the thin shell model (Fig. [T) at the 
ages < 6 X 10' yr. However, at the final stage of the ejecta decel- 
eration (f > 10* yr) the luminosity behavior differs from that of 
the thin shell model. Indeed at this phase the shocked gas cannot be 
treated as thin shell. After about 10* yr the reverse shock attains the 
center. This is accompanied by the overall contraction; as a result 
the emission measure increases and the luminosity attains maximal 
value ~ 3x 10'^ erg s"' . This is followed by the expansion which re- 
sults in the luminosity drop. At the most luminous phase, > 10'' 
erg s"', the SN la remnant lives ~ 4 x 10' yr. The temperature of 
the shocked ejecta at this phase is in the range of 10** - 10' K, so 
only a small fraction of the total luminosity (10 - 30%) falls into 
the standard Chandra band (0.2-10 keV). 

Generally, the SN expansion dynamics can be affected by the 
diffusion of relativistic protons into the ejecta. This process might 
modify dynamics by smoothing out the pressure jump at the bound- 
ary between ejecta and relativistic fluid. The time it takes to fill the 
SN by cosmic rays with the energy E per particle can be estimated 
as the time for the proton to escape from the relativistic bubble layer 
adjacent to SN. The volume comparable with SN of radius A is a 
spherical layer of a thickness of ~ Q.'iR. Assuming Bohm diffusion 
coefficient D = cTg/3 one gets the diffusion time 



(0.3Rf 
AD 



1.5 X lO'S, 



30 pc 



(15) 



The time it takes to fill the SN by cosmic rays at the essential de- 
celeration epoch turns out tremendous compared to the SN age 
(~ 3 X 10' yr). We conclude, therefore, that the diffusion pene- 
tration of relativistic protons into the SN envelope unlikely affects 
the SN expansion dynamics. 
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Time, yr 

Figure 4. X-ray luminosity of tlie reverse shocl^ in tlie model with the ICM 
temperature of Iff* keV (top) and the radii of the contact discontinuity and 
reverse shock (bottom). The luminosity decreases after ~ Iff* yr because of 
the significant expansion of the postshock layer between the reverse shock 
(bottom, lower curve) and the contact surface (upper curve). 

3.2 SN bulk motion 

After the expansion braking the SN shell still retains the bulk mo- 
tion with the typical velocity v, = 400 km s"'. If the deceleration 
of the bulk motion were negligible, the SN shell would escape the 
relativistic bubble after the average residence time Rtlv^ ~ 10^ yr. 
We now check, whether the ram pressure and the Alfven wave drag 
can substantially decelerate the bulk motion inside the relativistic 
bubble. 

The ram pressure drag force exerted on the SN shell is = 
nr\jPaV-, where the ambient density is pa = 3p/c-, assuming parti- 
cles dominate in the pressure. Using the equation of motion 

dv , , 

M— = -Txr^pav- , (16) 
at 

the characteristic deceleration length can be estimated as 

'"-5^-69^7''' kpc. (17) 

This shows that the deceleration of the bulk motion of SN by the 
ram pressure can be neglected. 

The Alfven wave drag exerted on the spherical SN shell 
is defined similarly to the case of the wind shell, i.e., Fj = 
(l/2)B^(v/vA)r^„- The chracteristic deceleration time is tj ~ 
v^MIFj, while the deceleration length, a v^t^j, is 

Id ~ — ~ O.SivooSj'pl''' kpc, (18) 

where Vj4oo = Vj/(400 km s"'). This shows that the Alfven wave 
drag eesentially brakes the bulk motion of the SN shell at the dis- 
tance much smaller than the radius of the relativistic bubble even 
for weak field 6 ~ 3 x 10"* G. We thus conclude that the ejecta of 
SN la exploded in the relativistic bubble cannot escape the bubble 
due to the bulk motion. Amazingly, the SN material is decelerated 



in the slow bulk motion at the distance ten times larger than in the 
maximum radius attained in the high speed envelope expansion. 



3.3 Rayleigh-Taylor instability and spike deceleration 

The swept-up SN shell decelerating in the light relativistic fluid 
is liable to the Rayleigh-Taylor instability (RTl) which generally 
should result in the fragmentation of SN shell, close to the stage of 
the significan deceleration, i.e., at about the stopping radius r^„. The 
situation is similar, albeit inverted with respect to Crab nebular. The 
Crab shell accelerated by the shocked relati vistic wind shows l ong 
thin RT spikes directed backward the center tester etal.ll996h . In 
case of decelerated SN dense RT spikes protruded forward could 
travel large distances before they stop. Yet it should be emphasised 
the difference with the Crab. In the latter case the SN material pres- 
surized by the relativistic plasma is cool with the thermal velocity 
"crab ~ 10 km s"' . The SN la material in the adiabatic reverse shock 
is hot with the thermal velocity u^^ ~ km s"'. The density con- 
trast in the Crab is therefore by factor (wsn/fcrab)^ ~ 10* larger. 

The behavior of RT spike may be affected by the KHI. For 
a cylinder spike of the radius a moving with the velocity v along 
its own axis the perturbation growth time for the most destructive 
wave number k ~ 1/a is fxH ~ (a/v);^^''^, where ^ is the density 
ratio of SN and bubble material. At the SN radius ^ = 20 pc 
the contrast ~ 10^. The distance at which the most dangerous 
mode grows is then ~ vty^ ~ ~ 300a. We do not aware of 

any multi-dimensional hydrodynamic simulations of a dense cloud 
moving in a rarefied relativistic fluid. A close analogue is the two- 
dimensional hydrodynamic simulations of a dense cloud moving 
in the post-shock intercloud rarefied gas l liaeinetal.ll 1994b . These 
simulations show that the cloud life time with respect to the frag- 
mentation and fragment deceleration is order ~ A(alv)x^^^, roughly 
four times larger than the Kelvin-Helmholtz time. Adopting this 
characteristic time one finds that the distance at which the spike will 
be destroyed and decelerated is ~ AvtKH ~ 4a^''^ ~ lO'fl. Assum- 
ing a /R ~ 10~^, co mparable with the fingers in case of Crab nebula 
(Hest er etal J 19961) . one finds that RT spike can travel ~ 10^ < 0.3 
kpc, rather small distanc e compared to the bubble radius. 

On the other hand, iNulsenI ^1982) argues that the increase of 
the width of the boundary layer due to the KHI quenches the very 
instability. As a result the mass loss is suppressed and eventually 
is defined by the momentum transfer [cf. equations ^ and Q]. In 
this case the problem of mass stripping due to KHI is reduced to 
the problem of a deceleration of RT spikes which is analysed in the 
next section. 

The longitudinal magnetic field also can suppress the KHI. 
Let pi and p2 be the density of rarefied and dense fluid. Accord- 
ing to the crite r ion of the KHI in the presence of magnetic field 
IChandrasekhail ( Il96 ih the condition that the magnetic field turns 
ofi^ the KHI is 

B>(27rpi)"\ = (67rp)"2(^j« 1.4x 10^* G, (19) 

where p = 10"'° erg cm"', c is the light speed, and v = 10' cm s"' 
are used. The required magnetic field B > 1.4 x 10"* G is within 
the range of field strength in the relativistic bubble, which can be 
as large as several 10"^ G. The magnetic stabilization of RT spikes 
against KHI thus seems quite plausible. Hereafter we address the 
deceleration of RT spike assuming its stability. 
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3.3.1 Drag in collisionless case 

A typical RT spike is assumed to be a cylinder with the mass m, 
radius a, and length b, which are assumed to remain constant. The 
RT spike pesumably moves along the axis in the relativistic plasma 
dominated by the relativistic particles. With the giroradius rg ~ 
3 X 10" cm and the RT spike radius a ~ lO^^r^ ~ 10"* cm only 
a motion along the regular magnetic field can be collisionless. The 
mean free path for relativistic protons propagating along the field 
can be constrained by scattering on the perpendicular component of 
a random field. Th e resulting mean fr ee path along the mean field B 
is ~ rg(BISBf dStrong et alj2007h . where 6B is the amplitude of 
a random field with resonance wave number k^^^rg ~ 1. Following 
IStrongetalJ l l2007h we assume the power law spectrum of random 
field energy density Wik) oc k^" with = 1.67, between maximal 
length scale k„,i„ ~ l/R/, and and minimal scale ki„a.x ~ 1 /'■.?• The 
integrated energy d ensity is normalized according to suggestion by 
IStrong et al] ( |2007|) : W = B-/(Sn). With these prerequisites one 
gets (6B/Bf ~ (k,„i„/k,,,y-^ ~ 10"^ and ~ lOV^ ~ 3 x lO'** 
cm, i.e., > Vg. The situation thus is about collisionless, although 
uncertainties of the relevant parameters do not preclude collisional 
regime as well. One needs therefore to consider both collisionless 
and collisional cases. 

For b » a the moment exchange occurs primarily via the cos- 
mic ray collisions with the lateral surface of the spike. The mo- 
mentum transferred to the colliding particle with the energy E as- 
suming diffusive reflection, is ~ (E/c^)v. For the particle flux on 
the unit of surface area (l/4)nc, where n is the cosmic ray concen- 
tration, the moment transferred per second to all the striking pro- 
tons (drag force) is Fd = (l/2)nabe(v/c), where e = 3p is the en- 
ergy density of relativistic particles. Note that the drag force could 
be derived also using average energy g ain of a rela tivistic parti- 
cle per collision with the cloud E(vlcf lFermi|[r949 h. Indeed, for 
a spherical cloud the total energy loss per second in that case is 
nc?ncE{vlc)^ = nc^ev^ jc. This implies the drag force nd^eivlc), 
which coincides with the above expression for Fj within the geo- 
metrical factor. 

The equation of motion of the spike in the collisionless case 
with the above value of Fj then reads 



dv V 
m— = -l.Snabp- 

dt c 



(20) 



We neglect here the head-on collisions which would contribute the 
term of the order ~ a/b <c 1 in the right hand side. Note, the 
transition from spike to the spherical blob corresponds to b = 2a 
in the drag force expression. The characteristic time of the spike 
deceleration is thus 



2 mc 
td~ - 



(21) 



3 nabp 

The ratio m/na- can be expressed via the surface density of the 
SN shell at the stopping radius as m/na^ = rjMIAnr]^, where the 
parameter » 1 because the spike is formed by a shell patch with 
the radius » a. The deceleration distance for the spike is then 



10 c a a 
Id = vtd ~ —Wsn-T = ^-^T]- kpc , 
9 V b b 



(22) 



where we make use of equation I ll2t and adopt v = 10' cm s"' 
and /■„, = 36 pc. For ?/ ~ 10 and ajb ~ 0. 1 the spike can travel 
~ 1 kpc before it gets completely decelerated. In collisionless case 
the RT fragments are decelerated efficiently inside the relativistic 
bubble. It should be stressed, however, that the collisionless regime 
can be realized only in the case of the motion along the magnetic 



field so, only a small fraction of RT spikes experiences this type of 
a deceleration. 



3.3.2 Drag in collisional case 

If the mean free path for cosmic ray protons is small /I <c a, one 
expects that the drag force should be proportional to v-. The gen- 
eral condition for that is the large Reynolds number Re> 10^. To 
estimate Re we adopt the spike radius a ~ 10"^ Tj,, ~ 10"* cm, the 
spike velocity v = 10' cm, and A = rg?& the mean free path for rel- 
ativistic protons. Assuming B = 10"^ G, i.e., r,, ~ 3 x 10" cm one 
gets Re= 3av/crg ~ 3x lO"". The condition Re> 10' thus is fulfilled 
for A <3 X 10*rg ~ lO""" cm, which is rather soft requirement. 

In the collisional approximation the drag force is Fj = 
3na-p(v/c)^. Following the recipe of the previous section one ob- 
tains the spike deceleration distance 



(23) 



For V = 10' cm s"' and rs„ = 36 pc one obtains Ij ~ llrj kpc, rather 
large value that exceeds bubble radius (5 kpc) even for modest value 
ofri ~ 1. The ram pressure drag thus almost does not decelerate RT 
spikes inside the relativistic bubble. 



3.3.3 Alfven wave drag 

The Alfven wave drag can operate for the large conducting body 
a > (v/c)rg ~ 10'" cm, which is easily met for RT spikes. Using 
the express i on for the power radiated in the form of Alfven waves 
iDrell etalJ ('l965')) one can write the Alfven drag force acting on 
the plasma spike with the radius a and length b as 



Fd = ab , 

27T Va 



(24) 



where va = 6/ ^j^npa is the Alfven velocity. Note, it is only the 
lateral surface area (» 2ab), that matters. Following arguments of 
the previous sections one gets the spike deceleration distance 



/7r\^/2 ac a , 
Id ^ 5rir,„ ^- j ^ - = S'/^Sj kpc 



(25) 



For fiducial model p = lO"'" erg cm"', B = 10"' G, v = 10' 
cm s"' one obtains Ij « brjiajb) kpc. This result shows that in 
case a = b the blob can travel the distance exceeding the radius 
of the relativistic bubble even for moderate values of 7/ > 2. For a 
long spike {b » a) deceleration is by factor bja stronger and the 
deceleration distance is accordingly shorter. For bla ~ 10 and 77 ~ 
10 the RT spike can travel ~ 6 kpc, a distance comparable with the 
adopted bubble radius Rh = 5 kpc. It takes roughly Rb/v ~ 5 x lO"" 
yr for the spike to reach the the bubble boundary assuming the spike 
average velocity of 10* km s"'. 

The effect of the Alfven wave drag is determined by the mag- 
netic field. If the field is weak, S = 3 x 10"* G, the deceleration 
distance is ~ 18 kpc, substantially larger than the bubble radius. 
On the other hand, for B > 10"' G the deceleration distance is 
smaller than the bubble radius and significan amount of RT frag- 
ments will remain in the relativistic bubble. The escape probability 
for RT spike in case of i? = 10"'' G can be estimated adopting mean 
free path A = vt^ = 6 kpc, i.e., t = Ri,/A = 5/6. Equation 1 11 It gives 
in this case the escape probability p^^c ~ 0.6. 
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Table 1. Parameters of fiducial model 



Parameter 


Description 


Value 


Rb 


Bubble radius 


5 kpc 


tb 


Age 


3 X 10' yr 


P 


Pressure 




B 


Magnetic field 


10-5 G 


n 


Number density of ICM 


0.02 cm-^ 


M, 


Stellar mass in the bubble 


9 X 10'" Mo 


N 


Stellar death rate 


0.83 yr-' 


Nsn 


SN la rate 


0.0066 yr-' 


rw 


Wind stopping radius 


0.5 pc 




SN stopping radius 


36 pc 



3.3.4 Spike deceleration by wind material 

The average number of wind shells along the average distance to 
the bubble boundary for the fiducial model is ~ 0.08 (Section[2j2} 
which means that for SN fragments the probability to collide with a 
wind shell is low, ~ 0.06. A question arises, what happens, anyway, 
if the collision takes place. 

The deceleration is determined by the ratio of column densi- 
ties (yu) of the projectile (RT spike) and target (wind shell). For the 
wind shell 

= = 1.4 X 10-' g cm-' , (26) 

where r„, = 1.5 x 10'** cm and m„ = 0.5 Mq are used. The column 
density of the spike ji^ = rjM/(4nr^^^) ~ 2 x lO^^rj g cm-- <c ;u„,. 
This comparison shows that the spike will be fully decelerated in 
a single collision with a wind shell. We conclude therefore that for 
the fiducial model there is non-negligible probability, ~ 0.06, that 
the spike will be decelerated in the bubble via collision with the 
wind shell. 



4 DISCUSSION 

The aim of this paper has been to get an idea on what happens to 
the matter ejected by the stellar population of BCG embedded into 
a bubble of relativistic plasma. We have found that the expansion of 
a wind envelope lost by a star is stopped by the pressure of the rela- 
tivistic fluid when the radius attains ~ 0.5 pc. Because of the small 
size of the wind shell its bulk motion is not decelerated neither by 
the ram pressure, nor by the Alfven wave drag, so the shell escapes 
the bubble together with the parent star. The SN la exploding in the 
relativistic fluid expands up to much larger radius ~ 30-40 pc and 
its bulk motion can be efficiently decelerated by the Alfven wave 
drag. Unless a RT instability operates, the SN material would not 
escape the relativistic bubble, but will instead be advected by the 
buoyantly rising relativistic fluid. 

The RT instability can strongly modify the behavior of the 
ejecta. In the framework of our fiducial model we find that signifi- 
cant fraction of SN fragments escapes the relativistic bubble. In our 
analysis of this scenario we rely on the fiducial model parameters 
outlined in Table 1. A gas lump crossing the bubble boundary and 
entering the ICM is decelerated after sweeping the ICM mass com- 
parable with its own mass. The escaping wind shell gets decelerated 
in the thermal plasma of ICM at the length of /„, ~ 1 kpc (Section 
I2.2I >. The wind deposits ~ 1.4 x 10' M© in this layer. This value 



should be compared with the mass of the ICM gas, which is already 
there. For R/, = 5 kpc and hydrogen number density n = 0.02 cm-^ 
the mass of the ICM in the boundary layer with the thickness /„, = 1 
kpc is 7 X 10' Mq, i.e., factor ~ 5 larger than the mass deposited 
by the wind shells. The deposited wind mass scales with the bub- 
ble radius as oc /J^, whereas the ICM mass in the layer is oc Rj^. We 
thus conclude that for the fiducial model the wind material escaping 
the relativistic bubble does not change significantly the density of 
the surrounding ICM; moreover the effect even smaller in case of 
M 87 in which bubble radius Ri; ~ 1 .4 kpc. The elfect of the energy 
deposition by the escaping wind shells is also negligible because 
the bulk velocities of the wind shells is subsonic and the deposited 
mass is low. The chemical composition of ICM is not affected by 
the escaping wind shells either. 

Unlike the wind shells, fragments of SN ejecta escaping the 
relativistic bubble may have a profound effect on the enrichment of 
ICM by iron peak elements. The resulting abundance is determined 
by the width of the mixing layer. Generally, one should consider 
time-dependent model of the formation of mixing layer. However, 
we assume that the mixing layer in the fiducial model is formed by 
the cumulative effect of all SNe. Emploing the momentum conser- 
vation arguments a deceleration distance for the RT spike entering 
the ICM with the velocity v, and decelerated down to v/ turns out 
to be 

M 

hn ~ J]Mvi/vf) ^ « 0.5j] pc , (27) 

snP 

where we used Vi/vf = 20, n = 0.02 cm-', M = 1.4 M©, and Tj,, = 
36 pc. For T] = 10 one gets the deceleration length 4n ~ 5 pc. At first 
glance the mixing layer could be identified with the deceleration 
layer. This layer contains M,„„ ~ 7 x 10^ Mq of the ICM gas. The 
mass produced by SN la during the life time of the bubble is M,„ ~ 
3 X 10^ Mq. If most of the SN mass escapes the bubble, the amount 
of escaping iron in the mixing layer turns out to be ~ 10' Mq which 
corresponds to the iron abundance ~ 80x (solar). 

On the other hand, the mixing layer could be broader because 
the total volume of shocked SN fragments in pressure equilibrium 
substantially exceeds the volume of the deceleration layer. Simple 
estimate based upon the pressure equilibrium suggests the total vol- 
ume occupied by shocked SN ejecta to be V, = Ns„E/p which im- 
plies the layer width ~ 300 pc. If this is identified with the width of 
the mixing layer than the iron abundance will be factor two larger 
compared to solar abundance of pre-existing ICM. The increase of 
the iron abundance by factor two changes the 0.6-2 keV emissiv- 
it>|by factors of 1.8, 1.5 and 1.3 for the gas temperatures 1, 2 and 
3 keV respectively. 

The latter estimates suggest complete mixing of injected iron 
with the ICM, which may not be the case. The point is that the de- 
celeration of SN fragments in the ICM results in the strong heating 
of the ejecta material up to temperature corresponding to its kinetic 
energy. Most of the iron therefore ends up in the high entropy/low 
density gas with the very little X-ray emission. The observational 
outcome thus critically depends on the mixing degree between hot 
SN gas and the relatively cool ICM. 

The above picture of mixing layer is very crude and one can- 
not rule out a possibility that the expanding relativistic bubble con- 
tinuously catches up with the mixing layer outer boundary so that 
most of the decelerated SN material injected in the ICM eventually 



0.6-2 keV is the energy range where present day grazing incidence X-ray 
telescopes are most sensitive 
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turns out engulfed by the relativistic fluid. Qualitatively this sce- 
nario then gets similar to the case when SN la ejecta bulk motion 
is decelerated well inside the bubble. In this case no strong X-ray 
emission is expected from fully expanded shells (because of the 
very low gas density). The fate of the iron generated by SN la in 
this scenario solely depends on the evolution of relativistic plasma, 
which can escape the central region of the galaxy , as suggested 
by ob servations of e.g. M87 jChurazov et alJuOOll ; iForman et al.l 
l2005l) . 



5 CONCLUSIONS 

We consider the outcome of the mass ejection by stars via winds 
and supernovae inside a bubble of relativistic plasma inflated by 
an AGN in the core of BCG. Wind shells are likely to escape the 
relativistic bubble and deposit their mass in the ICM within ~ 1 
kpc from the bubble boundary. SN la exploded inside the bubble 
is efiiciently decelerated owing to the pressure of the relativistic 
fluid. If the SN shells remain spherical until the expansion of the 
envelope stops and do not fragment, then they would not escape 
the bubble. In this case the iron produced by SN la is advected by 
the relativistic plasma and may leave the central region of the BCG 
together with buoyantly moving bubbles. 

As a possibility we consider a scenario in which the RT in- 
stability of the SN envelope at the deceleration phase breaks the 
shell into multitude of RT spikes. The analysis of the deceleration 
of RT spikes in the relativistic fluid shows that the SN fragments 
are able to escape the bubble. The fragments are decelerated in the 
ICM in a close vicinity of the bubble boundary thus producing Fe- 
rich layer. In the optimistic scenario this Fe-rich layer can enhance 
X-ray emission around bubbles of relativistic plasma, producing 
bright rims around bubbles. 
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